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LASER  PROPAGATION  TEMPORAL  SPECTRA 


I.  INTRODUCTION 

This  is  a  technical  report  under  Contract  F30602-72-C-0305 
titled  "Investigation  of  Laser  Propagation  Phenomena".  This  effort 
is  aimed  at  providing  theoretical  support  to  the  RADC  Laser  Propa¬ 
gation  Program  in  the  areas  of  atmospheric  propagation  phenomena  and 
microturbulence  statistics.  The  efforts  in  this  report  are  in 
support  of  the  experimental  program  being  conducted  at  RADC.  The 
immediate  aim  is  to  give  a  detailed  derivation  of  the  temporal 
spectrum  of  phase  difference  fluctuations. 

The  work  on  temporal  spectra  is  of  interest  because  it  provides 
information  on  time  scales  of  phase  and  arrival  angle  useful  to  systems 
designers,  and  further  phase  difference  fluctuations  are  intimately 
related  to  angle  of  arrival  fluctuations  for  small  apertures 
[Zintsmaster,  1971]. 

More  generally  temporal  spectral  measurements  provide  a  simple 
method  for  manipulating  experimental  data  to  give  a  check  on  our  knowledge 
of  the  light-beam  turbulence  interaction.  It  is  much  easier  to  instru¬ 
ment  time-spectral  measurements  because  data  need  be  taken  for  only  one 
apparatus  setting.  The  corresponding  spatial  correlations  or  spatial 
spectra  require  measurements  with  apparatus  adjusted  for  many  point 
separations . 

Temporal  correlations  and  spectra  are  in  themselves  useful 
quantities  to  compare  with  theoretical  estimates  to  gain  information 
about  atmospheric-turbulence  parameters.  In  this  context,  they  do  not 
suffer  from  aperture  limitations,  as  do  spatial  measurements,  because 
readily  available  time  lags  are  large  enough  to  represent  spatial 
separations  much  larger  than  those  conveniently  available.  Thus  they 
provide  a  convenient  optical  means  with  which  to  examine  outer-scale 
effects,  which  are  important  in  the  low  temporal  frequency  region  of 
the  phase  difference  spectrum.  Further  we  will  show  that  the  high 
temporal  frequency  dependence  of  the  phase  difference  spectrum  is 
easily  related  to  the  small  separation  dependence  of  the  phase  structure 
function,  thus  allowing  a  convenient  experimental  check  of  DS(P). 

The  calculation  will  be  done  in  three  steps.  In  the  first  step 
the  definition  of  time  lagged  phase  difference  spectra  is  written 
in  terms  of  the  phase  structure  function  using  Taylor's  frozen  tur¬ 
bulence  hypothesis.  The  result  is  in  terms  of  a  multidimensional 
integral  in  rectangular  coordinate  variables  whose  integrand  includes 
a  generalized  isotropic  Von  Karmen  index  spectrum.  The  second  step 
involves  evaluating  for  large  and  small  frequencies  all  of  the  indicated 
integral  ions  except  the  range  integration  to  obtain  the  differential 
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path  contribution.  The  result  is  put  in  terms  of  convenient  non- 
dimensional  variables.  The  third  step  is  to  perform  the  range  inte- 
q  rati  on  for  large  and  small  normalized  nondimensi  onal  temporal  frequencies 
and  thereby  obtain  limiting  forms  for  the  complete  solution.  These  forms 
are  tabulated  for  the  generalized  index  of  refraction  spectrum  and  for 
the  usual  k"11^  high  spatial  frequency  dependence.  The  complete 
solution  is  obtained  for  k inertial  subrange  dependence  by 
numerical  integration  for  several  normalized  separations  of  interest. 


II.  DERIVATION 


The  experimental  situation  is  shown  in  Fig.  1,  where  a  wind  of 


Fig.  1.  Diagram  of  the  experimental  measurement. 

constant  velocity  VT  is  blowing.  The  spherically  diverging  beam  is 
centered  on  the  z  axis.  The  +x  direction  is  defined  so  that 

(1)  VT  =  V  x  +  V„z, 

where  V  is  the  magnitude  of  the  transverse  wind  and  Vn i s  the  component 
of  the  wind  parallel  to  the  propagation  path..  Vn  may  be  either  a 
positive  or  a  negative  quantity.  The  phase  difference  (<J>2"‘f,l)  ^w0 

phase  points  separated  by  a  distance  p"  in  an^x-y  plane  is  recorded 
as  a  function  of  time.  The  vector  distance  p  makes  an  arbitrary, 
angle  e  with  the  +x  axis,  thus  enabling  us  to  determine  the  Fourier 
transform  (F)  of  the  time  laaged  phase  difference 

(2)  W^s(a),t)  =  F< ( 4>2 ( P2 * ) ”4‘i  (p-j  >i+t) )  (<J>2(p2’t)'<h  (pl  ,t^> 

for  all  relative  orientations  between  the  phase  points  and  the 
transverse  component  cf  the  wind  (Vx). 
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(3) 


( A-B )  (C-D)  =  ^(A-D)2-{A-C)2+(B-C)2-(B-D)2] 
in  Eq.  (2)  to  obtain 

(4)  W6$  =  \  f|<U2  ,r'-t )-^  (p,  ,t)]2-[*2 ,x+t)-^2 ,t)] 

+  [<J>1  (pi  , T+t )-4>2  (p"2  »'t)]2_C^'j  (pi  ,T+t)-(}li  (pi  ,t)]2j‘  • 

Each  term  in  Eq.  (4)  may  be  evaluated  in  terms  of  a  phase  structure 
function 


(5)  <  (<t>(p.j  >t.j ) ~4> (p  j >tj ))  >  -  (pi -p j  )  * 

A  general  expression  for  the  wave  structure  function  (Dw (p ) )  is 
obtained  by  multiplying  Hufnagel 's  (1964)  expression  for  the  square 
of  the  index  of  refraction  difference  by  the  wave  propagation  constant 
squared  (k2).  The  results  are 


(6) 


Dw^ 


=  2k* 


de 


n{Dn[0x  !’•> y 


where  £  is  the  difference  range  coordinate, n  is  the  average  range 
coordinate. 

The  phase  structure  function  varies  from  equality  to  one-half 
the  wave  structure  function  value  as  the  range  increases.  Hence 


(7) 


where  indicates  that  the  numerical  coefficient  is  between  0.5  and 

1.0.  Alternatively  Eqs.  (6)  and  (7)  are  derived  from  the  spectral 
expansion  solutions  for  the  wave  and  phase  structure  functions  in 
Appendix  A. 


Using  Taylor's  frozen  turbulence  hypothesis 

(8)  -  I^-Vyt,)). 
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Substituting  Eqs.  (5),  (6),  (7)  and  (8)  into  Eq.  (4)  yields 


F  n  'nf  k2  I  dn 


V1?x"Vt’CV5"Vzt)  +  Dn^Lpx“Vx» 


-  2Dn(-VT,0,c-V  x) 


assf!.  that  the  turbulence  is  stationary  so  that  Wfis(u)  is  not 
t.f™C}l0"0f .tme  on?ln  (t)-  Equation  (9)  may  be  written  using  the 


explicit  Fourier  transform  relationship 


«6s(“>  =  dx  Ds{x)  e'JU1T 


and  in  terms  of  the  index  of  refraction  spectrum 


00  00  00 

Dn(X,Y,Z)  =2  |  d<x  j  dKy  J  d<z  (l-cos(<xX+KyY+KzZ) 


-OO  -oo  -00 


x  Vw1^1 


The  result  is 


00  00  00  00 


WfiS^  =  2  \l!oj  k2  j  dx  |dn  j  dC  j  dlCx  {  dlCy  j  dlC2 


—  00  —  00  ..  00 


Cn(n)$IS0)(Kx’Ky’Kz) 


•j^exp  [j(<xRx+KyRy+KzC"E+T)}  'exp  {"J*(KxRx+KyRy+Kz5_e"'r)} 
-exp  |j(-KRx-KyRy+Kzc-E+T)|  -exp  {-j (-KRx-KyRy+Kzc-e_T)}- 
+2exp  |j(Kz5-e+T)|  +2exp  |-j(kz5-e"t)| 


L. _ A 


where  the  cosine  terms  have  been  put  in  terns  of  exponentials  and 
where 


R  h  pp  cose  is  the  range  dependent  x  component  of 
ray  separation 

R  prsine  is  the  range  dependent  y  component  of 
"  ray  separation 

£“  I^V+K^V  ±W  - 

Equation  (12)  is  the  complete  expression  for  the  time-delayed 
phase  difference  spectrum  and  hence  concludes  the  first  step. 

We  now  perform  each  of  the  indicated  integrations  except  that 
for  the  average  range  coordinate  n»  so  as  to  obtain  the  differential 
path  contribution  to  W(SS(oj).  The  £  integration  is  easily  performed 
resulting  in  six  delta  functions  involving  kz  in  the  remaining  integrand. 
The  k z  integration  is  performed  with  the  result 


(13) 


V“> 


■  {?-} 


where  g  =  k  V±w. 

X  X 


?rk  dx 


L 

/ 


00  00 


dn 


d<. 


d< 


-CO  -00 


-exp(j(K  R  +k  R  -e» 


Si^n  ^Kx,Ky’^|_  "r  £  v'x"x"'y  "y  “xwJ 
-exp|-j(KxRx+,cyRy-e'xT)}-exp[j(-KxRx-<yRy-^T) 

-exp{-j(-KxRx-KyRy-£;T)]+2exp{-jexT} 

r._i" 

+2exp'|+jexx  - 


Equation  (13)  shows  that  only  the  transverse  wind  velocity  (V)  con¬ 
tributes  to  the  temporal  spectrum,  which  agrees  with  Tatarski’s 
qualitative  analysis.  The  relatively  small  end  effects,  which  arrise 
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when  the  longitudinal  wind  (Vz)  moves  turbulence  into  and  out  of  the 
beam  at  the  transmitter  and  receive r>  are  insignificant  for  phase 
point  separation  much  less  than  the  path  length. 

The  above  result  is  a  direct  consequence  of  the  often  used 
approximation  of  extending  the  ^(=z2-zi )  integration  limits  to  ±»  in 
the  original  derivation  of  Ds(p). 

The  t  and  kx  integrations  are  done  in  a  similar  manner  with 
the  result  that 


(14) 


Wss(o)) 


AV1 


L 

► 

do 

0 


<u 


dK 

y 


—  tx> 


Cn(n)4n0>(ijf’Ky'0> 


where  <t>^^  is  even  in  each  of  its  arguments. 

Equation  (14)  may  be  put  into  a  convenient  form  by  defining 
dimensionless  variables 


(15a)  52  'f-  =  normalized  frequency, 


(15b)  R(n,v)  =  =  normalized  separation 

Lcrn' 


(15c) 


82(n,v) 


B(v)  V 
ujLq  (n ) 


1 


=  normalized  outer 
scale . 


(15d)  u  =  =  normalized  range. 


The  v  dependence  of  the  constant  B(v)  is  described  later  in  Eq.  (17c). 
In  terms  of  the  new  variables 
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06) 


6S 


2.2,  ..-1  -1 

it  k  LV  p  i} 


du 


Cn(u)3(u,v) 


f  <  Cl  -  ;jr  cos([-cose  +  Kssinejnu) 


0 

-  i-cos([cos0  +  >cssi ne]f«j)]  (-,^*,0) , 
c  n  p  p 


where  <  =  and  the  ever.sss  of  the  *  integration  has  been  used  to 

change  the  integration  range.  Equation  (16)  is  quite  general  because 
neither  the  functional  dependence  of  C^(u)  where  u  is  the  normalized 
average  range  coordinate  nor  the  form  of  (<x*Ky  *Kz)  has  been  specified. 

In  order  to  determine  the  differential  path  contribution,  we  now 
define  a  generalized  index  of  refraction  spectrin 


(17a) 


A(v) 


where  L0(n)  is  the  range  dependent  outer-scale,  and 


(17b) 


A(v) 


8tt 

0 (v+1  ) 


A(l  )  =  1/(2tt2) 


r(l-v)  si n[ (1  -V ) 


0<v<l 
1  <v<2 


(17c) 


B(v)  = 


■’(— ) 
,1-v  *  1  2  > 


r(Tf-) 


Equation  (17a)  is  a  generalized  Von  Karmen  index  spectrum  with  a  high 
spatial  frequency  dependence  of 
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2,  2  2 

K  +  ►  +► 

X  y  Z 


(v+3) 

2 


n 

T!,o  generally  used  value  of  v  is  -j  for  which  A^)  =  0.033005  and 
B(|)  =  1.0710. 

The  parameters  A(v)  and  B(v)  are  defined  so  that 


D„(r) 


n  o 


r  «  L, 


r  >>  L, 


Equations  (17b)  and  (17c)  are  derived  from  Eq.  (17a)  and  the  limiting 
forms  (Eq.  (18))  in  Appendix  B. 

We  now  substitute  Eq.  (17a)  into  Eq.  (16)  and  after  rearranging 
terms  obtain 


r~  r(^) 
VJL  £ 

2  r(^) 


(v+2)  [  du  C2(u)[f(v+2)(u,v) 


X  [  dx[l-  lcos([-cose+Ke(u,v)sine]fiu) 


-  goos([cose+K3(u,v)sine]nu)] 


X  (1+K  ) 
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When  | (±coso+3sinG)fiuj  <<  1  Eq.  (19)  becomes 

,  W.s(a)(i)  [>r(^-r 

(20) _ °S  P  .ll  VTT  2  I  _-v 
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0.5 

1.0 


i  5/2 (r (l2)).,  .  2,  v+1 


L2  r(^)J 


duC^(u)e"(v+^)(u,\))^— 


.2  „ ,  2  2/  x  .  2 


c1k[cos  0-t-K  3  (u,v)sin  e][l+K  ] 


2-,  1  2  ' 


The  k  integration  is  easily  performed  giving 


(21a) 


w6s  («)£)_ 


0.5 


r/yi2\ 

5/2/ n  2  M  * /  %,  2,  v+1 

’  l^a)A(',)l<  Lp 


=  ft 


-v 


1  _ 
f 


du 

o  ■- 


cos2e+g  ^u,v^  sin2e 

V 


<(u)s(u,v)-(“+2)(^). 


and  for  a  horizontal  path 


(21b) 


W5s<<) 


0-515/2/  r(  2  Ar(  x.2.  v+1  r2 

1-0J  i^i+3 -JA(^)k  LP  Cn 


(o2+r2)  (  2  } 


X  [n2cos2e+(^L)sin2e]  fdu  ~ 

V  J  0 

0 


and 

(21c) 


8ir5/2  ffll76)  A(2/3)  =  4 -384. 


9 


When  | (±cosG+3sine)ou|  <<  1  Eq.  (19)  becomes 

-1 


(20) 


0.5 

1.0 


l„V2[<l]A(v)kV+l 

rm/ 


~  r  /v+2  v-| 

/IT  r(~' 


L2  r(^)J 


-V 


0 


duc2(u)S-(v+2)(u,v)i 


-(41) 


dK[cos29+K2s2(u,v)sin2e][l+K2]  2 


The  k  integration  is  easily  performed  giving 


(21a) 


0.51  S/2/r(  2 1  |  ..2.  v+1 

l.or  l7^±^JA(v)k  Lp 


=  sf  v 


du 


cos2e+- sin2e 

V 


d(u)3(u,v)“(v+2)(^-). 


and  for  a  horizontal  path 


(21b) 


=  (o2+r2)  ‘  2  1 


Lpv+1C2 
n 


and 

(21c) 


X  [a2cos2e+(5^-)sin29] 


du 


3u 


A(2/3)  =  4.384. 
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Thus,  there  is  a  quadratic  variation  (along  with  any  or 
e(u,v)  variation)  of  the  differential  path  contribution  to  W6s(u)  with 
range  when  Q  is  small. 

When  | (cose±3(u,v)sinG)fiu|  »  1  the  cosine  terms  oscillate  rapidly 
and  we  replace  them  by  their  average  value  (0).  The  phase  difference 
speccrum  is 


(22a) 


W,.(a)£) 

6S  p 


v+2 

A(v)k2Lpv+^ 


-  n-(v+2) 


1 


=  n 


8 


fo.5 

1.0 


duC^(u)s‘^v+2^ (u,v) 


0 


1  r(V+3) 

V'  2  ’> 


and  for  a  horizontal  path 


(22b) 


8i 


0.5 L  5/2 

li.or 


IvCf-Y 
lr<#> 


!A(v)k2Lpv+1C2 


-(— ) 
(n2+R2j  2 


du 


Equations  (22)  show  that  for  large  un,  the  only  variation  in  the  dif¬ 
ferential  path  contribution  is  due  to  c£(u)  or  g(u,v).  Figure  2  shows 
the  limiting  forms  and  the  break  point  location  when  the  path  is 
horizontal.'  A  typical  differential  path  contribution  curve  would 
smoothly  pass  from  the  low  frequency  form  to  the  high  frequency  form 
and  then  oscillate  with  rapidly  decreasing  amplitude  about  the  high 
frequency  (constant)  amplitude.  This  completes  the  second  step. 

The  third  step  is  to  integrate  Eqs.  (21b)  and  (22b)  to  obtain 
the  final  equations.  The  resultant  equations  are  tabulated  in  Table 
1,  for  the  generalized  index  spectrum  and  in  more  detail  in  Tables  2, 

3  and  4  for  the  usual  case  of  v=2/3.  In  each  of  the  four  tables  the 
spectral  amplitude  has  been  multiplied  by  2  and  the  range  of  u  halved  to 
give  0<to<°°(one  sided  spectrum).  In  order  to  get  the  complete  horizontal 
path  solution  for  W5S(u),  we  change  the  order  of  integration  in  Eq. 

(19)  and,  assuming  constant  and  L0,  perform  the  u  integration  to 
obtain 
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TABLE  I 


The  limiting  forms  of  the  one  sided  phase  difference  spectrum 
for  any  relative  orientation  (e)  of  phase  points  and  wind 
velocity  (V),  and  for  a  Von  Karmen  refractive  index  spectrum 
having  an  index  structure  function  Dn  %  Cfirv  when  r  <<  L0. 


B(v) 


1/v 

0<v<2 


2,  v+1 


l(n2+R2)"(v+2)/2  n2cos2e+  ^-sin2 


±ncos 


+(n2+R2 )^2sinel  >> 


.  (o2+R2)-^+2)/2 
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TABLE  II 


The  one  sided  phase  difference  spectrum  limiting  forms 
for  any  orientation  of  phase  points  and  wind  velocity 


n  = 


2  Trfp 


n  .  1 . 071 P 
R  =  — j- 

0 


| tftco s  +(n2+R2 )^2sine  1  << 


u  =  2/3 


z„2 . 02 v-4/3  [Vcos20  A  (a2+R2)sin2o 
>/3  “  R  '  L  6  4 


8-7H?’oPl  c-p5/3 


_  r“8/3  H  COS  0  +  R  sijl .6.  s  R  >>  f/ 

6  4 


9  9 

„-2/3  cos^e  j  sinS~| 

n  l_“T~  4  J  ' 


n  >>  R 


l+ncos  +(fl2 +R2)1/2sine|  >>  /S 


=  (n2+R2)“4/3 
-  R'8'3 


R  »  n 
fi  »  R 
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TABLE  III 


Phase  difference  spectrum  limiting  form  when  the  wino 
direction  is  parallel  to  the  line  connecting  the  phase 
points 


=  hh. 


R  = 


LOTk 

L 

o 


6  =  0‘ 


R  <<  /G 


v  =  2/3 


R  >>  n 

SG  >>  n  >>  R 
n  >>  & 


R  >>  /G 


,  1  r-8/3  n2 


G 


n 


-  »-8/3 


637>V77  »  a 
a  »  6371V77 
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/ 

y 


tabu:  iv 


Phase  difference  spectrum  limitinq  forms  when  the  wind 
direction  is  perpendicular  to  the  line  connecting  the 
phase  point 


R  — 7-^  ,  o  =  90°  ,  v  =  2/3 


R  <<  /6" 


W  ,(n)t-) 

6  S 


;^k2L  c2p5/3 


1  n-2/3 
4  K 


=  1  o'273 


R  M 

2  -«-*  il  R 


q  »  2 


R  >> 

•  R-8/3 


R  »  n 
n  >>  R 
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J 


(S:)(-v-) 

OS  L' 


r  ('’?■) 


-JO.51 1  ''2"''  _5/2.,  ,k2.  v+lr2 
8] ^  ri  ,o  ^  A(v)k  Lp  C 


i;(x> 

2  ,  /  v + 3  \ 

•  \~o~i 


(v+2) 


x  d< 


1_  sinR-cose+K/Ssineifl. 
2  {-  cose+KSsineln 


l  -  (— -) 

1  sinRcose+KBsinoln]  r-,+  2-,  2 

?  { cos Q+xBsin Gift  L  J 


The  integrand  in  Eg.  (23)  contains  the  spatial  frequency  filter 
functions  and  the  normalized  index  of  refraction  spectrum.  Figure  3 
shows  the  filter  functions  for  a  particular  normalized  separation  (R) 
and  normalized  frequency  (u),  and  three  orientation  angles  (o).  The 
normalized  index  spectrum  is  also  presented.  For  0°  and  90'  the  two 
sinx/x  terms  combine  to  yield  a  single  term.  The  O'  case  has  a 
constant  amplitude  dependent  only  on  u.  The  901  case  filter  function 
is  a  sine  function,  which  discriminates  against  the  low  index  of 
refraction  spatial  frequencies.  The  general  o  case  filter  function 
is  represented  by  e=15°  which  has  two  1/2  amplitude  sine  functions 
that  are  displaced  from  the  origin. 

The  integral  in  Eq.  (23)  is  readily  found  using  numerical 
integration.  The  zero  degree  case  is  easily  evaluated  because  the 
filter  factor  is  a  constant.  In  this  particular  case 


V<> 


jo.5).  im , 

Vo;  re£> 


5/2  A(,)k2Lpu+1C2 


-(“ ) 

5^][R2+ft2]  2 

SI  J 


Figures  4-7  present  the  computer  generated  curves  for  R  =  10,  1.0, 
0.1,  0.01;  for  o  =  0°,  5°,  15°,  45°,  90°  and  v=2/3.  Also  included  as 
narrower  lines  in  the  figures  are  the  asymptotic  limits  determined  from 
Table  2  (e  =  0°,  90°  can  also  be  obtained  from  Tables  3  and  4). 
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FILT 


ft  ~  NORMALIZED  FREQUENCY 


Normalized  phase  difference  spectrum  for  a  horizontal  path 
versus  normalized  frequency  and  orientation  angle,  when 
the  normalized  separation  is  1.0.  The  light  lines  are  the 
limiting  forms  from  tables  2,  3  and  4. 


Fig.  6.  Normalized  phase  difference  spectrum  for  a  horizontal  path 
versus  normalized  frequency  and  orientation  angle,  when 
the  normalized  separation  is  0.1.  The  light  lines  are  the 
limiting  forms  from  tables  2,  3  and  4. 
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10  1  10  •  1.0 

NORMALIZED  FREQUENCY 


Normalized  phase  difference  spectrum  for  a  horizontal  path 
versus  normalized  frequency  and  orientation  angle,  when 
the  normalized  separation  is  0.01.  The  light  lines  are  the 
limiting  forms  from  tables  2,  3  and  4. 


Figure  4  is  for  the  case 


R  =  =  10 

o 


which  represents  a  large  phase  point  separation  compared  to  the  outer 
scale.  The  low  frequency  offset  between  the  computed  curves  and  low 
frequency  asymptote  is  due  to  a  higher  order  term  (a  constant)  that  is 
not  included  in  the  low  frequency  asymptotic  solution.  An  interesting 
feature  of  this  figure  and  following  Figs.  5,  6  and  7  is  that  the 
normalized  spectrum  approaches  a  constant  level  for  small  normalized 
frequencies  except  for  the  case  when  the  wind  is  exactly  parallel  to 
the  phase  point  separation  (e=0°).  A  second  general  result  is  that  the 
nigh  frequency  spectra  is  independent  of  orientation  angle  (0)  and  the 
normalized  separation  (R),  and  hence  it  is  also  independent  of  the 
outer-scale . 


It  is  the  mid  frequency  range 
normalized  separation.  For  small 


R< 


e  =  0°  1 
e  =  90°  p 


that  is 


characterized  by  the 


the  mid  frequencies  which  have  a  -2/3  slope  extend  from  ft=R  to  n-/6 
for  0=0°  and  to  ft=2  for  e=90°.  Thus  for  small  R,  that  is  phase  point 
separations  much  less  than  the  outer  scale,  the  lower  mid  frequency 
break  point  determines  R  and  hence  L0  since  the  phase  point  separation 
is  known. 


We  may  also  determine  the  exponent  v  in  the  phase  structure 
function  (Eq.  (18))  from  phase  difference  spectrum  data  by  taking  the 
ratio  of  parallel  to  perpendicular  spectra  amplitudes  for  the  case  of 
small  normalized  separation  ;  (R)  and  small  or  intermediate  nonnalized 
frequencies  (ft).  From  the  equations  for  the  limiting  forms  given  in 
Table  1,  the  ratio  is 


This  ratio  is  accurate  when  ft  <<  R  <<  1  or  R  <<  ft  <<:  1  as  can  be  seen 
in  Fig.  7  which  shows  both  the  limiting  form  results  and  the  exact 
results  calculated  for  the  particular  value  of  v  =  2/3.  The  in¬ 
equality  (R«ft<l )  defines  the  more  useful  ft  range  because  the  ratio 


22 


simplifies  to  itself.  Therefore  we  have  a  convenient  means  of 
testing  the  fundamental  index  of  refraction  structure  function 
equation  which  is  widely  used  in  turbulence  theory  (see  Eg.  (18)). 


HI.  PHYSICAL  PICTURE 

In  order  to  get  a  better  physical  picture  of  the  factors  in- 
fluencing  the  phase  difference  temporal  spectrum  we  re-examine  several 
of  the  previously  developed  equations. 

First  put  the  fundamental  definition  of  the  phase  difference 
spectrum  given  in  Eq.  (2)  in  terms  of  correlation  functions, 

(25)  w,ss  ('■’)  =  '(2Bs(0,r)  -  B$  (♦„,•,)  -  B^-,,,,)). 

Hence  W*s  („,)  is  the  Fourier  transform  of  a  combination  of  four  corre¬ 
lation  functions.  The  first  two  correlations  which  are  combined  into 
one  term  are  between  the  phase  of  one  ray  (see  Fig.  1)  and  the  phase 
of  the  same  ray  at  a  time  x  later.  This  is  equivalent  to  the  same  time 
correlation  between  two  parallel  rays  separated  by  a  distance  vx  arrisinq 

cress  ^  °r  %fr°2en  3hurbul6nCe  TOdel-  "he  last  terns  are  phase  9 

rav  Lh  Ih!  h  "h'i  Th^  Cr°SS  corre1atlons  are  between  the  phase  of  one 
y  and  the  time  delayed  phase  of  the  other  ray.  Again  usinq  Tavlor's 

frozen  turbulence  hypothesis  the  time  delay  correlation  “  equivalent  to 

tran^vPrspdliPlaH-men!-°f  the  UpWard  ray  in  a  diction  opposite  to  the 
s  a  -p1  d  nIn  one  case  the  rays  are  moved  further  apart 

?°  ?aL  uhey  do  not  cross  (°i2lL)  and  in  the  other  case  they  approach 
each  other  or  cross  at  a  point  intermediate  to  the  ray  end  points. 

Now  consider  the  differential  path  contribution  to  Ws.(„)  from 

Eq"  n ?21h Th  ?LE?-  (2b)’  as  refl6Cted  on  the  corresponding  terms 
(21b;  a"d(22L)-  The  flrst  ten"  is  Proportional  to  range  and 
gives  a  constant  differential  path  contribution  while  the  last  terms 
have  a  pronounced  range  dependent  variation  due  to  the  ranqe  dependent 

hrc  ;°tribu  tW°  rayS-  tHe  ran9e  te™  in  the  differential 

path  contribution  to  W6s  increase  with  range  and  tend  to  approach  a 

constant  toward  the  end  of  the  path,  thus  qualitatively  verifying  the 
range  dependence  of  Eqs.  (21b)  and  (22b).  y  y  9  th 

refrarti vp^inHp v°cin+^Sxi9ate+the  physica1  interaction  between  the 
sider  Jhp%  nf  s?atlal  spectrum  and  the  experiment  geometry  con- 
^dPr  thk •  ^"dimensional  phase  grating  model  of  the  atmosphere  where 

qratinamwi?htl0n-°J-Va1heS  °f  and  <z  “responds  to  a  phase 

greting  with  periodic  phase  variation  in  the  <  direction.  The 

_z  clirected  gratings  influence  each  ray  in  exactly  the  same  wav  so 

<7  In  Fan°min^  UtS  Phase  difference  and  hence  the  zero  for 

z  hq.  (]9j  (for  example).  Hence  we  consider  the  qratinqs  which 
have  normals  in  the  x-y  plane.  Figure  8  further  shows9  Examples 


MEASUREMENT 

POINTS 


(a  ) 


(b) 


Phase  grating  model  of  the  experimental  situation 


of  grating  orientations  (with  <.7=0),  that  is,  gratings  with  s-  directed 
perpendicular  and  parallel  to  trie  separations..  The  gratings  with 
normals  perpendicular  to  the  separation  do  not  influence  the  phase 
difference  while  the  parallel  case  has  the  maximum  effect  on  the  phase 
difference.  Hence  the  filter  function  of  Eq.  (19)  has  projected  the 
gratings  with  normal  in  x-y  plane  onto  the  separation  vector. 

It  is  the  x  component  of  ttie  index  spectrum  which  couples 
spatial  and  the  temporal  frequency  variation  of  the  phase  difference 
as  the  wind  moves  past  the  phase  points.  Hence  it  is  the  >x  =  u/V 
component  that  produces  the  frequency  w  in  the  temporal  spectrum. 
Furthermore,  if  a  particular 


[-->] 

(n  an  integer)  which  has  an  associated  gratong  wavelength,  n>x,  that  is 
equal  to  the  local  x  component  of  the  ray  separation  (Rx),  then  the  two 
rays  experience  identical  phase  fluctuations  and  thus  zero  phase  difference 
at  that  point  on  the  ray.  Hence  there  is  a  zero  in  the  differential  path 
contribution  to  the  temporal  spectrum  at  that  point.  Maxima  on  the 
temporal  spectrum  occur  when  n  is  replaced  by  (n+1/2).  This  effect  is 
found  in  Eq.  (14),  upon  expanding  the  cosine  of  the  sum  of  two  angles  and 
combining  terms.  These  maxima  and  minima  are  in  the  plane  wave  spectra 
because  there  is  only  one  separation  for  all  points  on  the  rays.  But 
the  maxima  or  minima  at  one  point  on  the  spherical  wave  path  is  added 
to  the  contribution  at  other  points  where  the  filter  function  is 
neither  at  a  maxima  or  a  minimum  and  hence  the  spherical  wave  spectra 
is  nearly  free  of  local  oscillations  in  amplitude. 

As  a  final  comparison  we  note  that  the  low  temporal  frequency 
spectrum  approaches  a  constant  amplitude  because  of  u/V  =  kx  has  selected 
a  spatial  frequency  in  the  saturation  region  of  the  index  spectrum.  The 
exception  is  the  0°  case  where  the  filter  function  has  a  strong  influence. 
The  high  frequency  temporal  spectra  goes  as 


(-) 

V 


-(v+2) 


which  results  from  projecting  the  ^  index  spectrum  in  the  p 

direction.  This  completes  the  description  of  salient  features  of  the 
temporal  spectrum  in  terms  of  physical  concepts. 
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IV.  CONCLUSIONS 


The  phase  difference  temporal  spectrum  has  been  evaluated  for 
arbitrary  orientation  of  the  phase  observation  points  in  a  plane 
perpendicul ar  to  the  beam  axis.  The  calculation  has  been  based  on 

pnfiflH nnc dh°f  sPectral  expansions  and  the  Rytov  approximation  but  the 
2  !  haf  been. ^arranged  to  emphasize  their  relationship  to  the 
physical  picture  using  a  phase  screen  model.  The  calculation  uses  a 
generalized  index  of  refraction  spectrum  which  has  arbitrary  power 
la*  dependence  ,n  the  inertial  subrange  and  includes  outer-sca°e 
ertects.  This  generalization  allows  one  to  directly  test  phase 
difference  spectra  measurements  to  check  the  generally  accented  *-11/3 

outer-scale  ^  development  also  allows  a  direct  evaluation  of  the 


For  a  (<)  to  the  (v+3)  power  law  inertial  subrange  spatial 
frequency  dependence,  a  (<)  to  the  (v+2)  power  law  high  frequency 
phase  difference  temporal  spectrum  is  predicted,  independent  of 
phase  observation  point  orientation.  The  low  frequency  temporal 
spectra  are  strongly  orientation  dependent.  These  features  are 
indicated  in  graphs  of  typical  cases. 
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APPENDIX  A 


In  this  appendix  we  obtain  the  phase  structure  function  in  terms 
of  spatial  coordinate  integrations  from  the  often  used  spectral  ex¬ 
pansion  results. 

We  start  the  calculation  by  comparing  the  spherical  wave  phase 
structure  function  (Carlson,  1969) 


(Ala) 


and  the  spherical  wave,  wave  structure  function 


(Alb) 


8k27T2 


where 


p 

L 

cn<"> 


and  <  =. 


is  the  source  wave  number 
is  the  measurement  separation 
is  the  propagation  range 

is  the  turbulence  strength 


<Jl^(k)  is  the  turbulence  spectrum 


2  .  2 
K  +  K 

x  y 


The  phase  structure  function  in  Eq.  (Ala)  is  identical  to  the 
wave  structure  function  in  Eq.  (Alb)  except  that  the  former  has  an 
additional  cosine  squared  term  in  its  integrand.  For  the  case  of 
small  separations  (p)  and  long  ranges  (L)  the  cosine  squared  term 
oscillates  rapidly  during  the  significant  part  of  the  <  integration 


UhPnWe  Jdn1replace  ttle.cosine  squared  term  by  its  average  value  (0  5) 

When  p  is  large  and  L  is  small,  the  cosine  squared  term  remains  essentially 

?S2?  WAsni'SK  r;,°f  ss.;  :sn  SKS 

s  s.1  contrlbution  of  the  “,sr  s" 


(A2a) 


D  (p). 


(A2b) 


Using  Eqs.  (Alb)  and  (A2a),  and  setting  <’  =  <  —  yields 

n 

L 


Ds(p)  =  8k2,2  |“'5|  j  dr|  c2(n) 

0 


dK‘  <' 


f  subsequent  equations*  drop  the  prime.  <  is  the  maanituHp 

dependence  of  the^sna?^!1?  ^  X_Y  E13"6*  We  n**in  the  Mentation 
oepenaence  of  the  spatial  frequency  by  using  the  identity 

n 


(A3) 


27  |  <**(l-cos(^cos*))  -  l-Jo(USfi.) 


where  <j>  is  the  angle  between  <  and  p  in 
Eq.  (A3)  into  Eq.  (A2)  and  changing  to 


the  X-Y  plane.  Substituting 
rectangular  coordinates  yields 


(A4) 


L 

dn  C2(n) 

1  I 


00 

r 


-00 


(<X>Ky .0) 


^Kxpx+Kypy 


The  relationship  between  the  three  dimensional 
spectrum  and  the  two  dimensional  spectrum  is 


index  of  refraction 
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APPENDIX  B 


We  will  now  derive  the  constants  in  a  generalized  isotropic 
Von  Karmen  index  spectrum  such  that  the  corresponding  index  structure 
function  is 


(Bl) 


D„(rt  •  •{ 


r2  v 

Cn  r 


C2  Lv 
n  o 


r  << 


r  »  Lq 


where  C2  is  the  turbulence  strength,  and  L0  is  the  outer  scale. 

The  inclusion  of  the  parameter  v  allows  the  experimentalist  to 
test  the  widely  used  2/3  power  (v=2/3)  variation  of  with  separation  (r) 

A  general  expression  for  the  three  dimensional  isotropic  index  of 
refraction  spectrum  is 


(B2) 


v+3 

2 


where  the  constants  A(v)  and  B(v)  are  to  be  evaluated, 
has  been  chosen  such  that 


(B3) 


D„(r)  -  2 


1  I  1 


(l-cos(K-r)) 


*n(<) 


di< 


The  exponent 


The  general  inverse  transform  relationship  (Tatarski  1961)  for  iso¬ 
tropic  turbulence  is 


(B4) 


Dn(r)  =  8n 


00 


f 

,  sin  <r  i 

l 

0 

k  r  j 

Substitute  Eq.  (B2)  into  Eq.  (B4)  to  obtain 


d<  . 


(B5)  Dn(r)  =  8*  A{v)  Cj; 

. 


sin  <r\ 
k  r  j 


2 

K 


2 

+  K 


v+3 
2  2 


Then  change  variables 


(B6)  oj  =  <r 

and  obtain 


(B7 )  Dn(r)  =  8tt  A(v )  Cj;  rv 

. 


sinoo 

0) 


,  v+3 

Brr.  21  2 

g  “ 


When  r/L0  <<  1  we  drop  the  Br/L0  term  as  being  small  in  comparison  to 
w  during  the  significant  part  of  the  integration,  and  usinq  Eq.  (B1 ) 
obtain 


(B8) 


A(v)  = 


8tt 

L  0 


<-+1  > 

w  j 


that 


The  integration  is  put  into  a  more  tractable  form  by  noting 


(B9) 

. 


sinu  -(v+1  )  , 

-  a)  doo  = 

0) 


00 

■ 

a 

dw 

da 

c 

db  sin(buj)  u'v 


Equation  (B9)  is  easily  evaluated  when  v=l,  the  result  being  tt/4. 
When  v?n,  interchange  integration  limits  and  using 
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oo 


(BIO) 


sin(w)  u_v  h.j  =  r(l -v)sin  j(l-v) 


to  obtain  A(v)  when  0  <  | v- 1  |  <  1.  The  final  result  for  A(v)  is 


-1 


(B1!)  A(v)  = 


vf^+TT  r(1_v)sin(1“v'  2 


0  <  v  <  1 

1  .<  v  <  2 

v  =  1 


2tt 


When  U  =  k/B  is  substituted  into  Eq.  (5)  the  result  is 


(B12) 


Dn(r)  =  8ttC^A(v) 


★ 


u2du _ _ 

(1+u2)^3T/2 


When  rB/L0  >>  1  we  may  neglect  the  sine  term  and  easily  integrate  to 
obtain  an  expression  involving  the  remaining  constant  B(v).  Using 
Eq.  (Bl) 

(B13)  B(v)  =  [4tt  B ( 3/2  ,  v/2)  A(v)]v  . 


Upon  substituting  Eq.  (Bll  )  into  Eq.  (B13)  we  obtain  the  final 
expression  for  B(v). 


(B14) 


B(v) 


.(1-v) 


r(^) 


r(*£) 


1 

V 


The  complete  expression  for  the  phase  structure  function 
is 
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where  we  have  used  Tatarski's  (1961)  transform  pair. 


